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ABSTRACT 



Whether gamma-ray bursts are highly beamed or not is a very important 
question, since it has been pointed out that the beaming will lead to a sharp 
break in the afterglow light curves during the ultra-relativistic phase, with 
the breaking point determined by T ~ 1/9q, where T is the bulk Lorentz fac- 
tor and 9 is the initial half opening angle of the ejecta, and such a break is 
claimed to be present in the light curves of some GRBs. In this paper we will 
examine whether all the observed breaks in GRB afterglow light curves can 
be explained by jet effects. Here we present a detailed calculation of the jet 
evolution and emission, and have obtained a simple formula of bulk Lorentz 
factor evolution. We show that the light curves are very smoothly steepened 
by jet effect, and the shape of the light curve is determined by only one param- 
eter — (S/n) 1 / 8 ^^ 4 , where E and n are the fireball energy and surrounding 
medium density. We find that for GRB990123 and GRB991216, the jet model 
can approximately fit their light curves, and the values of (-Ssi/ru) 1 / 8 ^ 4 
are about 0.17 and 0.22 respectively While for GRB990510, GRB000301c, 
GRB000926 and GRB010222, their light curves cannot be fitted by the jet 
model, which suggests that the breaks may be caused by some other reasons, 
jet should be not the unique reason. 

Key words: gamma rays: bursts — shock waves — ISM: jets and outflows 



© 0000 RAS, MNRAS 000, 000-000 



Can all breaks be explained by jet effects? 3 

1 INTRODUCTION 

Multiwavelength follow-up of gamma-ray burst afterglows has revolutionized GRB astron- 
omy in recent years, yielding a wealth of information about the nature of GRBs (Klose 2000; 
Castro-Tirado et al. 1999 and references therein). The afterglows can well be explained as 
the emission from a relativistic blast wave which decelerates when sweeping up interstellar 
medium. The dynamical evolution of GRB fireballs and the emission features have been 
studied by many authors (e.g. Sari 1997; Meszaros, Rees & Wijers 1998; Wei & Lu 1998a, 
2001a; Huang et al. 2000a; Sari, Piran & Narayan 1998; Wijers, Rees & Meszaros 1997). 

One of the most important questions to the nature of GRBs is the total energy released 
in the event, which depends on two factors: the distance scale and the opening angle of GRB 
ejecta. The distance scale of several GRBs have been derived through measuring the redshift 
in optical observations, while the question whether the gamma-ray emission is isotropic or 
not remains uncertain. The observations show that, if the gamma ray emission is isotropic, 
then for some GRBs their total energy will be too large, for example, the isotropic energies 
of GRB971214 and GRB990123 are 3 x 10 53 ergs and 3.4 x 10 54 ergs respectively (Kulkarni et 
al. 1998, 1999). Such a crisis encountered by this extreme large energy forced some people to 
think that the GRB emission must be highly collimated in order to reduce the total energy. 

Then how can we tell a jet from an isotropic fireball? Rhoads (1997, 1999) has pointed out 
that the lateral expansion of the relativistic jet will produce a sharp break in the afterglow 
light curves, and such breaks are also claimed to be present in the light curves of GRB990123 
(Kulkarni et al. 1999; Castro-Tirado et al. 1999), GRB 990510 (Harrison et al. 1999; Stanek 
et al. 1999), GRB000301c (Rhoads &Fruchter 2001; Masetti et al. 2000) and GRB000926 
(Sagar et al. 2001; Price et al. 2001), GRB991216 (Halpern et al. 2000), GRB010222 (Masetti 
et al. 2001; Stanek et al. 2001). Here we give a detailed analysis of the jet evolution and 
emission under the relativistic case to examine whether the jet model can account for all the 
observed breaks in GRB afterglows, and find that for four bursts, GRB990510, GRB010222, 
GRB000926 and GRB000301C, their afterglow light curves cannot be explained by jet effects. 
In next section we discuss the jet evolution and emission in the relativistic regime, in section 3 
we fit six GRBs' afterglow light curves, and finally we give some discussions and conclusions. 
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2 THE EVOLUTION AND EMISSION OF JET 

Here we follow our previous paper (Wei & Lu 2000b) to calculate the emission flux from 
the jet. We assume that the radiation is isotropic in the comoving frame of the ejecta, the 
radiation cone is uniquely defined by the angular spherical coordinates (9, 0), here 9 is the 
angle between the line of sight (along z-axis) and the symmetry axis, and is the azimuthal 
angle. Because of cylindrical symmetry, we can assume that the symmetry axis of the cone 
is in the y — z plane. In order to see more clearly, let us establish an auxiliary coordinate 
system (x', y', z') with the z'-axis along the symmetry axis of the cone and the x' parallel the 
x-axis. Then the position within the cone is specified by its angular spherical coordinates 
9' and 0' (0 < 9' < 9j, < (f)' < 2n, 9j is the jet opening half angle, which increases with 
time). It can be shown that the angle between a direction (9', <p') within the cone, and 
the line of sight satisfies cos© = cos9cos9' — sin9sin9'siri(p'. Then the observed flux is 

F v = #' / sin9'd9'DZl'{vD- 1 )— (1) 
Jo Jo a 2 

where D = [T(l - pcosQ)}' 1 is the Doppler factor, (3 = (1 - T" 2 ) 1 / 2 , v = ZV, I'(u') is the 
specific intensity of synchrotron radiation at u', and d is the distance of the burst source. 
Here the quantities with prime are measured in the comoving frame. For simplicity we have 
ignored the relative time delay of radiation from different parts of the cone. 

It is widely believed that the electrons have been accelerated by the shock to a power 
law distribution n e {j) oc 7~ p for 7 min < 7 < ^max, and consider the synchrotron radiation 
of these electrons, we can obtain the observed flux 

F v oc f(T)zT 3 n^ +5 ^ +1)/4 u-^ 2 (2) 

where T is the time measured in the observer frame, f(T) = r^^ p+7 ^ 2 (T—lY P+1 ^ A y(TY P+5 ^ A ^f~^f3 3 , 
y(T) = 7 is just the ratio of specific heats, ^ min = £ e (T - l)^;^Ef, m p (m e ) is the 
mass of proton (electron), and when the line of sight is along with the jet symmetry axis, 
z = /5 _1 [(1 —f3)-(p +9 y 2 — (1 — f3fij)~ < - p+9 ^ 2 ], /ij = cos9j, n is the surrounding medium density, 
£b is the energy fraction occupied by magnetic field. In the relativistic case, we have 
F u oc T 2p+6 T 3 [1 - ( L^)(P+9)/2] n (p+5)/4 e (p+l)/4 (3) 

1 + r 9j 

Now let us discuss the jet evolution. For an adiabatic relativistic jet expanding in sur- 
rounding medium, the evolution equation of energy conservation is 

T 2 Nm p c 2 = E (4) 
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Figure 1. The evolution of FBq with T*. The circles are the numerical results, and the solid line is our fit results (eq.(6)). 

For constant density, the total particle numbers N = nV = 27rnr 3 (l — cos9j)/3 ~ nnr 3 9j/3 
for 9 j ; <§C 1. In the relativistic case, r ~ AT 2 cT, so we have 

r = 372(-^)V8^7 1 /4 T -3/8 (5) 

where E 5 i is the fireball energy in units of 10 51 ergs, ni is the surrounding medium density 
in units of 1 atom cm" 3 . The jet half opening angle 9j = 9 + 9' ~ 9 + f^ff , where 
c s is the expanding velocity of ejecta material in its comoving frame, and for relativistic 
expanding material it is appropriate to take c s to be the sound speed c s = c/3 1 / 2 (Rhoads 
1997, 1999). Then we can obtain the evolution of jet Lorentz factor T9q(1 + f^f f^) 1 ^ 4 = 

372(^i) 1 /8^/4 T -3/8_ go we gee that? j f the yalue of (^1)1/8^3/4 . g ^ ^ calculate 

the variation of T9o with time T, and here we obtained a simple fitted formula 

T9 = 5.2T- a(T * } (6) 

where T* = (^)~ 1 ^ 3 9Q 2 T day: T day is the observed time in units of 1 day, and a(T+) = 
0.375 + 0.075 il^/gooo^oL ■ It should be noted that the index a is a function of T*, so the 
relation between T9 and T* is not simply determined by a, the actual decay index k should 
be calculated as k = d\n{Y9o) / dhiT* = — (a + In T+da jd In T*). Fig.l gives the numerical 
results of T9q evolution and our fitted results. Equation (6) can be easily used to evaluate 
the evolution of bulk Lorentz factor, and we can derive the value of (E/n) 1 ^ 8 9^ 4: through 
fitting the afterglow light curve. 
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Figure 2. The R-band light curve of the GRB990123 afterglow, the contribution of host galaxy has been subtracted. The solid 
line represents our best fit. The fit parameters are: (E 51 /n 1 ) 1 / 8 B^ A = 0.17 ± 0.08, p = 2.49 ± 0.25, the reduced \ 2 = 2 -2. 

3 FITTING RESULTS 

Based on the results described above, now we fit six GRBs' afterglow light curves in 
which the breaks are present, they are GRB990123, GRB990510, GRB000301c, GRB000926, 
GRB991216 and GRB010222, the results are as follows. 

3.1 GRB990123 

GRB990123 was the brightest GRB seen by BeppoSAX to date, it is in the top 0.3% of all 
bursts if ranked by the observed fluence, its redshift is estimated to be z = 1.6, corresponding 
to a luminosity distance fi L ~ 12 Gpc, and the isotropic 7-ray energy w 3.4 x 10 54 ergs 
(Kulkarni et al. 1999). The observed spectral slope between the optical band and the X- 
ray wavelengths was fiox = —0.68 ± 0.05, corresponding to the electron distribution index 
p between 2.3 and 2.5 (Castro-Tirado et al. 1999). Fig.2 illustrates our best fit to the R 
band light curve, the contribution from the host galaxy has been subtracted. The best fit 
parameters are: {E^/n^e^ = 0.17 ± 0.08, p = 2.49 ± 0.25, and the reduced \ 2 = 2.2. 
We see that jet expansion can explain the observed light curve. 

3.2 GRB000926 

GRB000926 was a burst lasted about 25 s, its redshift is about z = 2.066, yielding the 
luminosity distance di = 16.9 Gpc, and the isotropic energy E n ~ 2.5 x 10 53 ergs. The value 
of spectral index in the X-ray -optical region is about (3 = —0.8 or -0.9, corresponding to 
the electron distribution index p between 2.6 and 2.8 (Sagar et al. 2001). Fig. 3 gives our best 
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Log t (days since burst) 

Figure 3. The R-band light curve of the GRB000926 afterglow, the contribution of host galaxy has been subtracted. The solid 
line represents our best fit. The fit parameters are: (-Esi/ni) 1 / 8 ^' 4 = 0.087 ± 0.053, p = 2.6f ± 0.36, and the reduced \ 2 = 6. 

fit to the R band light curve, the contribution from the host galaxy has been subtracted. 
The fitted parameters are: (.E^i/nx) 1 / 8 ^ 4 = 0.087±0.053, p = 2.61 ±0.36, and the reduced 
X 2 = 6, so we see that the jet model cannot fit the light curve well. 

3.3 GRB991216 

GRB991216 was one of the brightest 7-ray bursts detected by BATSE, with a fluence of 
2.56 x 10~ 4 ergs cm -2 above 20.6 KeV (Kippen 1999). Its duration is about 15 seconds, and 
the redshift is about z = 1.02, yielding the isotropic gamma-ray energy about 6.7 x 10 53 
ergs. The observed spectral slope between the optical band and the X-ray wavelengths was 
Pox = —0.74 ±0.05, corresponding to the electron distribution index p between 2.4 and 2.6 
(Halpern et al. 2000). Fig. 4 illustrates our best fit to the R band light curve, the contribution 
from the host galaxy has been subtracted. The best fit parameters are: (E 5 i/ni) 1 ^ s 6^ 4: = 
0.22±0.04, p = 2.51±0.08, and the reduced x 2 — 1-5. We see that jet expansion can explain 
the observed light curve well. 

3.4 GRB990510 

GRB990510 was also a strong burst with duration of about 100 s, its redshift is z = 1.619 ± 
0.002, corresponding to the luminosity distance di — 12 Gpc, and the isotropic energy 
E 1 ~ 2.9 x 10 53 ergs (Vreeswijk et al. 1999). The observed optical spectral index was f3 = 
—0.61 ±0.12, corresponding to the electron distribution index p between 2.1 and 2.3 (Stanek 
et al. 1999). However, we find this burst cannot be fitted by simple jet model, since the 
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Figure 4. The R-band light curve of the GRB991216 afterglow, the contribution of host galaxy has been subtracted. The solid 
line represents our best fit. The fit parameters are: (.Esi/rti) 1 / 8 ^'' 4 = 0.22 ± 0.04, p = 2.51 ± 0.08, and the reduced x 2 = 1-5. 
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Figure 5. The R-band light curve of the GRB990510 afterglow. The solid line represents our best fit. The fit parameters are: 
(E 5 i/n 1 ) 1 / H el ,i = 0.096 ± 0.004, p = 1.87 ± 0.03, t c = 438 ± 25 days, the reduced X 2 = 7.8. 

reduced y 2 = 23.5. One possible way to improve the fitting results is that the cooling 
frequency u c crosses the observed band around break time. Similar to eq.(3), we can obtain 
the observed flux for v > v c in the relativistic case 

^ocr 2 ^T 2 [l-( rT i^)^ 8 )/ 2 ] (7) 

Using eqs.(3)(7) and T9 evolution, we re-fit the afterglow light curve of GRB990510, the 
parameters are: {E^/n^e^ = 0.096 ± 0.004, p = 1.87 ± 0.03, t c = 438 ± 25 days (where 
t c is the time T* when the cooling frequency v c crosses the observed band), the reduced 
X 2 = 7.8. From Fig. 5 we see that the jet model still cannot explain the light curve well. 
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Figure 6. The R-band light curve of the GRB010222 afterglow. The solid line represents our best fit. The fit parameters are: 
(-Esi/ni) 1 / 8 ^ 4 = 0.169 ± 0.018, p = 2.01 ± 0.04, t c = 105 ± 35 days, the reduced x 2 = 5.1. 

3.5 GRB010222 

GRB010222 was detected simultaneously by the Gamma-Ray Burst Monitor (GRBM) and 
the Wide Field Camera (WFC) instruments on BeppoSAX, the 40 - 700 KeV fluence was as 
bright as (9.3 ±0.3) x 10~ 5 ergs cm -2 , only be surpassed by GRB990123. The redshift of this 
burst was determined at z = 1.477 (Jha et al. 2001). The spectral index of optical data is (3 = 
—0.88 ± 0.10, corresponding to the electron distribution index p between 2.5 and 3 (Stanek 
et al. 2001). However, this burst also cannot be fitted by simple jet model, since the reduced 
X 2 = 15.7. If we consider the case that the cooling frequency v c crosses the observed band 
around break time, we obtain the following fit parameters: (.E^i/ni) 1 / 8 ^ 4 = 0.169 ± 0.018, 
p = 2.01 ± 0.04, t c = 105 ± 35 days, the reduced \ 2 = 5-1- From Fig.6 we see that the 
jet model can approximately fit the light curve. However, it should be noted that the fit 
parameter p is obviously smaller than that from spectra measure, and in addition, Masetti 
et al. (2001) has pointed out that the constancy of the optical spectral shape excluded that 
the v c has crossed the optical band during the observations, so we think that the jet model 
still cannot explain the light curve. 

3.6 GRB000301c 

GRB000301c was a peculiar burst with duration of about 10 s, there are short term flux 
variability in its optical light curve, and besides this, the overall light curve shows sharp break 
(Sagar et al. 2000; Masetti et al. 2000). Its redshift is z = 2.0335 ± 0.0003, corresponding to 
the luminosity distance = 16.6 Gpc, and the isotropic energy E 1 ~ 3.4 x 10 53 ergs. The 
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spectral index is about (3 = —0.8 or -0.9, corresponding to the electron distribution index 
p between 2.6 and 2.8 (Sagar et al. 2000). However, this burst cannot be fitted by simple 
jet model, since the reduced x 2 — 31. If we consider the cooling frequency v c crosses the 
observed band around break time, the reduced x 2 value is still very high, x 2 — 25, so this 
burst also cannot be explained by jet model. 

4 DISCUSSION AND CONCLUSION 

Whether gamma ray bursts are beamed or not is a very important question, since it is related 
with the energy sources of GRBs, and the afterglow observations provide a very good chance 
to study this question. Rhoads (1997, 1999) pointed out that the sideways expansion of the 
relativistic jet would cause the blast wave to decelerate more quickly, leading to a sharp 
break in the afterglow light curve. However, some numerical calculations show that such 
break is much weaker and smoother than the prediction (Moderski et al. 2000; Panaitescu 
& Meszaros 1999; Kumar & Panaitescu 2000). We have also pointed out that the steepening 
of the afterglow light curves can be observed only when the beaming angle is very small, 
#o < 0.1 (Wei & Lu 2000b), but in that paper we fixed the Lorentz factor r = 300. 

In this paper we reanalyse the jet evolution and emission in detail, and we find that the 
effect of jet expansion on the afterglow light curve can be described by only one parameter 
— (E/nY^O^ 4 . From equations (3) and (7) it is obvious that, when T6j 3> 1, the observed 
flux F v oc T- 3 ^ 1 )/ 4 for v < v c or F v oc T^" 2 )/ 4 for v > u c , while when Y6j « 1, the flux 
F v oc T" 3p / 4 for v < v c or F v oc T~^ 3p+1 ^ 4 for u > v c if no lateral expansion, and F v oc T~ p if 
lateral expansion is important, these are the same as the previous analytical results, which 
claim that there will be two breaks in the afterglow light curve. However, from equations 
(3) and (7) we see that the afterglow light curve is steepened gradually and smoothly, which 
depends on the variation of T9j. 

The breaks predicted by theoretical models have been observed in some GRBs' afterglow 
light curves, and have been generally considered as evidence for collimation of the relativistic 
GRB ejecta. We have fitted six GRBs' afterglow light curves in which the breaks are present, 
and find that, for GRB990123 and GRB991216, their light curves can be approximated fitted 
by the jet model, but for GRB990510, GRB000301c, GRB000926 and GRB010222, their light 
curves cannot be fitted by the jet model. So we conclude that, although the lateral expansion 
of the relativistic jet can lead to a break in the afterglow light curve, this may be not the 
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unique reason, there should be some other reasons to cause the break. For example, the 
transition from relativistic phase to non-relativistic phase of the blast wave may cause such 
break (Dai & Lu 1999; Huang et al. 2000b); the effects of inverse Compton scattering can 
flatten or steepen the light curves (Wei & Lu 1998b, 2000a); the curved emission spectra 
(not a simple power-law as usually assumed) can also produce a break in the light curves 
(Wei & Lu 2002). 

From fitting we have obtained the values of (-E^i/rii) 1 / 8 ^ 4 . As an approximation, this 

3 /4 

value can also be regarded as the value of 6 , since it is very weakly dependent on the value 
of (£51/711). Therefore from the fitting we can estimate the opening half angle of GRBs, for 
GRB990 1 23 6 « 0.094, and for GRB991216 O « 0.13. In addition, from equation (5) we 
see that, if the value of (E^/riiY^Ol^ is known, then we can also get the value of T 6 t^ 8 , 
here t is the deceleration time. In a few bursts, soft X-ray emission has been observed from 
the end of the GRB phase, indicating that the external shock had already set in by the end 
of the GRB (at i 7 ) (Giblin et al. 1999). In other cases there is no detectable X-ray emission 
after the GRB, suggesting that t 7 < to (Pian et al. 2001). In order to constrain T 9o, we 
assume that the observed GRB duration is a good measure of t (see also Panaitescu & 
Kumar 2001), and equation (5) shows that T has a moderate dependence on t , thus the 
error due to this assumption is likely not too large. Furthermore, using the values of 9 
derived above, we can give an estimate of the initial Lorentz factor, we find for GRB990123 
To ~ 130, and for GRB991216 To ~ 228. These values of Tq seem reasonable since it is 
widely believed that the initial fireball Lorentz factor should be larger than 100 in order to 
avoid photon-photon annihilation. 

Here we only consider the case that the line of sight is along the jet axis, since only this 
can give an analytical results. For the case that the line of sight is off the jet axis but still 
within the cone, we have calculated the light curves using equation (1), and found that the 
overall shape of the light curve is similar to the case 6 = 0, but the break is more smooth, 
as suggested by Ghisellini & Lazzati (2000), who did not consider the lateral expansion. 

In summary, our conclusion is that, not all the breaks in the afterglow light curves can 
be explained by jet model, so the steepening of the light curves may be caused by varied 
reasons, jet should be not the unique reason. But if the light curve can be fitted by jet model, 
then we can obtain the value of (E 5 i /ni) 1 / 8 ^ 4 , and furthermore we can give a constraint 
on the values of r and 9q. 
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